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A  Numerical  and  Experimental  Investigation  of  Sabot  Separation  Dynamics 


Martin  J.  Guillot,  William  G.  Reinecke 


1.0  INTRODUCTION 

High  acceleration  loads  transmitted  to  gun  launched  projectiles  during  the  launch  cycle 
require  projectiles  to  be  structurally  supported  while  in-bore.  This  is  even  more  critical  for  long 
rod  penetrators  under  investigation  today.  The  usual  method  for  supporting  projectiles  in-bore  is 
to  use  a  sabot  made  of  two  or  more  petals  which  encases  the  projectile  and  provides  the  needed 
structural  support  for  the  projectile  to  withstand  the  acceleration  loads.  However,  after  the  projec¬ 
tile  h^  exited  from  the  muzzle  of  the  gun,  the  sabot  must  be  removed  from  the  projectile  as 
quickly  as  possible  to  minimize  aerodynamic  drag.  This  is  normally  accomplished  by  designing 
the  sabot  so  that  aerodynamic  forces  strip  the  sabot  petals  from  the  projectile.  Sabot  separation  is 
a  highly  complex  and  unsteady  aerodynamic  process,  and  if  not  accomplished  properly  can  be  a 
source  of  perturbing  forces  and  moments  that  can  result  in  an  unacceptable  level  of  dispersion 
Therefore,  a  thorough  understanding  of  sabot  separation  dynamics  has  been,  and  continues  to  be' 
an  important  area  of  research  for  the  U.S.  Army. 

In  an  effort  to  develop  the  capability  to  predict  sabot  separation  dynamics  within  the 
framework  of  a  reasonably  simple  physical  and  computational  model  capable  of  rapid  parametric 
predictions,  the  AVCO  Systems  Division,  under  contract  to  the  U.S.  Army  Ballistics  Research 
^boratory  began  developing  a  computer  model  in  the  late  1970’s  to  use  as  a  sabot  design  tool. 
This  code  has  become  known  as  the  AVCO  Sabot  Separation  Code.  It  was  originally  written  to 
compute  sabot  separation  for  relatively  simplified  geometries.  However,  the  innovative  projectile 
designs  and  increasingly  complicated  sabot  geometries  under  investigation  today  have  limited  the 
usefulness  of  the  current  version  of  the  AVCO  code  when  modeling  these  sabots.  Additionally,  the 
onginal  AVCO  code  has  undergone  several  modifications  over  the  years  by  different  reseL-ch 
agencies.  These  modifications  have  been  documented  in  a  series  of  reports  [1-6].  However,  to  the 
authors’  knowledge,  there  is  no  single  report  that  documents  the  theoretical  development' of  the 
current  version  of  the  AVCO  code  entirely,  nor  provides  a  comprehensive  user’s  manual.  There¬ 
fore,  the  current  effort  was  undertaken  for  two  reasons.  The  first  was  to  modify  the  original 
authors  version  of  the  code  to  compute  trajectories  of  more  complicated  sabot  geometries  under 
development  today.  The  second  was  to  provide  comprehensive  documentation  and  a  user’s  man¬ 
ual  for  the  current  version  of  the  code. 

.u  T  modified  and  original  versions  of  the  AVCO  code  were  run  for  sabots  launched  at 

thclAT  Hypervelocity  Launch  Facility.  The  results  of  both  versions  of  the  code  were  compared  to 
X-Ray  and  witness  plate  data  of  radial  and  angular  displacement  at  specific  downrange  locations 
to  assess  the  overall  ability  of  the  code  to  accurately  predict  the  sabot  trajectories,  and  to  assess 
the  improvements  of  the  modifications  to  the  original  version  of  the  code. 

In  addition  to  the  modifications  to  the  AVCO  code,  a  relatively  simple  analytical  model 
was  developed  to  help  verify  the  trends  predicted  by  the  AVCO  code.  The  results  of  the  analytical 
model  are  also  compared  to  the  experimental  data  and  are  used  to  identify  possible  further 
improvements  to  the  aerodynamic  model  in  the  AVCO  code. 
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2.0  THEORETICAL  AND  NUMERICAL  FORMULATION  OF  AVCO  CODE 
2.1  Equations  of  Motion 

The  theoretical  formulation  of  the  model  incorporated  into  the  AVCO  code  was  o  la: 
developed  by  Crimi  and  Siegelman  [1].  In  keeping  with  the  intention  to  provide  a  single  ap 
hensive  documentation  of  the  AV CO  code,  and  to  provide  the  reader  with  a  base  refere’  ^  point 
for  the  discussion  of  the  modifications  made  in  this  study,  the  original  formulation  b  isented 
here. 

The  basic  equations  of  motion  governing  the  sabot  petal  and  projectile  dynam  s  are  the 
conservation  of  linear  momentum  (Newton’s  second  law)  and  the  conservation  of  angular 
momentum.  These  equations  are  written  in  an  inertial  coordinate  system,  (X,Y,Z),  for  each  seg¬ 
ment  and,  for  asymmetric  discard,  the  projectile.  The  inertial  coordinate  system  is  chosen  such 
that  the  origin  is  at  the  center  of  gravity  of  the  projectile,  which  is  assumed  to  move  at  constant 
velocity.  The  centerline  of  the  projectile  lies  along  the  X-axis  with  the  positive  direction  toward 
the  projectile  tip.  See  Figure  1.  Also  shown  in  the  figure  are  body  fitted  coordinates  used  for 
computing  the  aerodynamic  forces.  Initially,  the  X  and  z  axes  are  aligned  with  each  other.  The 
equations  of  motion  are  written  in  vector  form  as: 


YF  =  mX 


(2) 


(Hg);c>-z  change  of  angular  momentum  measured  in  body  fitted  coordinates, 

(x,y,z),  CO  is  the  angular  velocity  of  the  sabot  segment  measured  in  (X,Y,Z)  coordinates,  is  the 

angular  momentum  of  the  segment,  and  the  vector,  X,  represents  the  three  cartesian  coordinate 
components. 

Equations  (1)  and  (2)  are  sufficient  to  describe  the  motion  of  the  sabot  petals  and  projec¬ 
tile.  However,  it  is  convenient  to  define  additional  variables,  and  to  solve  the  equations  in  terms  of 
those  variables.  First,  the  rectangular  components  of  the  angular  velocity,  co ,  are  related  to  the 
time  rate  of  changes  of  Euler  angles,  (0,  (j),  \(/)  ,  which  are  also  shown  in  Figure  1.  Referring  to 
the  figure,  it  is  easily  seen  that  the  rectangular  components  of  angular  velocity  are  related  to  the 
time  rate  of  change  of  the  Euler  angles  by  the  matrix  relation. 


cos(t)  sin0sin<t)  0 

0 

-sintj)  sin0cos(t)  0 

V 

= 

0  COS0  1 

A 

Solving  (3)  for  the  Euler  angles  yields. 
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Figure  1.  Sabot  and  inertial  coordinate  system. 
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o)^sin(|)  -  (o^cos<t) 
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(4) 


9  =  co^cos(j)-co^sin<j) 

=  ci)^- cot  0  (cousin  cos  (j)) 

Second,  the  components  of  angular  momentum  with  respect  to  the  sabot  petal 
given  by  the  following  matrix  relation. 


(5) 

(6) 

center  of  gravity  are 
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(7) 


Equation  (7)  is  used  to  write  the  rectangular  components  of  the  angular  velocity,  m ,  in  terms  of 
the  components  of  angular  momentum,  {h^,  h  h^) ,  in  equation  (2)  and  equations’ (4)  through 
(6).  This  yields  a  system  of  coupled  equations  given  by  equations  ( 1)  and  (2),  and  (4)  through  (6) 

These  equations  are  solved  using  the  Runge-Kutta  method,  which  is  described  in  more  detail  in 
Section  3.0. 

The  aerodynamic  forces  are  most  easily  evaluated  in  body  fitted  coordinates,  but  before 
they  can  be  used  in  the  above  equations  they  must  be  transformed  into  inertial  coordinates.  The 
transformation  is  given  by. 
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Py 

where,  the  transformation  matrix  [  T]  is  given  by, 

(cos(t)cos\(t- sin(t)cos9sin(j))  (- sin(l)cos\|/- cos(t)cos9sin\|/)  (sin' 

(cos(|)sin(|)  + sin())cos9cos\i/)  (- sin<l)sin\|/+ sin(l)cos9cos\|/)  (-sin9cos(t)) 
(sin9sin(t))  (sin9cos(t))  (cos  9) 

2.2  Aerodynamic  Model 

The  original  aerodynamic  model  developed  for  the  AVCO  code  consists  of  two  primary 
flow  regimes.  The  first  is  when  the  sabot  petals  are  in  close  proximity  to  one  another  and  the 
forces  and  moments  on  each  segment  are  influenced  by  the  presence  of  the  other  petals,  and  by 
shock  reflections  and  interactions.  This  is  termed  the  interaction  region.  The  second  regime  com¬ 
mences  when  each  sabot  petal  is  in  free  flight  and  not  affected  by  the  presence  of  the  other  sabot 
petals.  In  the  interaction  region,  a  relatively  simple  aerodynamic  model  based  on  inviscid  2-D  gas 
dynamic  relations  is  applied  along  the  centerline  of  each  sabot  petal  to  determine  the  pressure  dis¬ 
tribution.  Newtonian  pressure  distribution  is  assumed  to  exist  when  the  sabot  petals  are  in  the  free 
flight  regime.  Forces  are  computed  using  both  the  interaction  region  procedure  and  the  Newtonian 
method  at  every  point  in  the  trajectory.  Whenever  the  forces  computed  by  the  interaction  method 
fall  below  those  computed  by  the  Newtonian  method,  the  sabot  petal  is  deemed  to  be  in  free  flight 
and  only  Newtonian  forces  are  computed  thereafter.  The  sabot  geometry  modeled  in  the  original 
AVCO  code  consists  of  a  conical  front  scoop  and  a  cylindrical  body.  The  relevant  geometric 
parameters  are  shown  in  Figure  2. 


SIDE  VIEW  FRONT  VIEW 

Figure  2.  Original  sabot  geometry  model  [5]. 
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In  the  interaction  region  the  pressure  distribution  is  divided  into  two  main  areas  of  the 
sabot:  the  sabot  front  scoop  and  the  sabot  underside.  On  the  front  scoop,  the  procedure  used  to 
compute  the  pressure  distribution  depends  on  the  angle  of  attack  of  the  sabot.  If  the  turning  angle 
of  the  flow  encountering  the  front  scoop  is  small  enough  for  an  attached  shock  to  exist,  the  pres¬ 
sure  on  the  front  scoop  is  assumed  to  be  constant  at  the  value  behind  the  oblique  shock.  Planar 
shock  relations  are  used.  If  the  turning  angle  is  greater  than  for  an  attached  shock,  the  pressure  on 
the  front  scoop  is  assumed  to  vary  parabolically  from  stagnation  at  the  outer  radius  to  sonic  at  the 
inner  radius.  The  two  situations  are  sketched  in  Figures  3  (a)  and  (b).  The  base  pressure  level  on 
the  underside  of  the  each  sabot  petal  is  computed  by  assuming  that  the  flow  turns  through  a 
Prandtl-Meyer  expansion  of  angle  (ti  -  0 ^).  However,  the  pressure  along  the  underside  of  the 

sabot  increases  in  the  axial  direction  from  the  tip  to  the  tail  due  to  shock  reflections  and  shock 
interactions.  An  approximate  model  based  on  empirical  relations  and  experimental  data  by 
Schmidt  [7]  is  used  to  approximate  the  pressure  distribution.  The  experimental  data  indicate  that  a 
pressure  pulse  occurs  on  the  underside  of  the  sabot  due  to  the  interaction  of  the  projectile  bow 
shock  and  the  reflection  of  the  sabot  bow  shock  off  of  the  projectile.  This  pressure  pulse  is  approx¬ 
imated  by  the  procedure  shown  in  Figure  4.  A  line  from  the  sabot  front  scoop  is  extended  parallel 
to  the  front  scoop  down  to  a  point  where  it  intersects  the  projectile.  A  Mach  wave  is  then  reflected 
back  up  to  the  sabot  underside  at  the  Mach  wave  angle,  |i .  This  serves  to  define  the  beginning  of 
the  pressure  pulse,  5^ ,  and  the  end  of  the  pressure  pulse,  ,  as  shown  in  the  figure.  The  expres¬ 
sions  for  5^  and  are  given  by. 


St 


dcosQj 

sin(0y.-a) 

S  I  ^sin0^  rcosacotp.- sina~[ 
*  sin  (0^-a)  Lsinacot|i+ sinaj 


(10) 


The  pulse  is  assumed  to  be  triangular  with  a  peak  halfway  between  5^  and  .  The  magni¬ 
tude  of  the  pressure  pulse  is  determined  from  turbulent  separation  data. 
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Figure  4.  Location  of  pressure  pulse  on  sabot  underside. 

The  validity  of  this  approach  is  discussed  thoroughly  by  Crimi  &  Siegelman  [2].  The 
authors  compare  the  pressure  distributions  predicted  by  this  model  with  experimental  data.  They 
achieve  good  agreement  between  their  predictions  and  the  experimental  data.  Based  on  this,  they 

conclude  that  this  approach,  while  having  little  theoretical  basis,  represents  a  valid  empirical 
model.  ^ 


3.0  MODIFICATIONS 

In  the  current  work,  several  modifications  were  made  to  the  original  formulation  presented 
in  Section  2.0.  Briefly  summarizing,  they  are: 

•  The  Runge-Kutta  algorithm  was  modified  to  allow  any  order  Runge-Kutta  scheme  to  be  used. 
It  is  only  necessary  to  input  a  table  of  coefficients  through  a  user  subroutine. 

The  error  estimator  and  time-stepping  scheme  were  improved  so  that  the  new  time  step  is 
computed  based  on  the  error  calculated. 

The  Newtonian  formulation  was  changed  from  straight  Newtonian  to  modified  Newtonian. 

The  same  side-face  areas  used  to  compute  the  forces  and  moments  in  the  interaction  region  are 
now  used  to  compute  the  forces  and  moments  in  the  Newtonian  regime. 

•  The  front  scoop  geometry  and  aerodynamic  model  were  modified  so  that  the  front  scoop 
geometry  can  be  modeled  with  one,  two,  or  three  line  segments. 

Each  of  these  modifications  is  discussed  in  more  detail  below. 
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3.1  Runge-Kutta  &  Time  Stepping  Scheme 

The  equations  presented  in  Section  2.1  are  solved  using  the  Runge-Kutta  r  lod.  e 
equations  are  rewritten  as  a  system  of  twelve  first  order  differential  equations  by  defi-'  i  a  vet  >r 

X’y  =  [X,Y,  Z,  6,  \j/,  (|),  9,  \j/,  (j),  h^,  hy,  h.]  where  7=1, 12 .  This  yields  a  sy<^  n  of  twelve 
equations  for  the  twelve  unknown  variables  which  are  written  as: 


=  ^4 


X2  =  Xg 


X,  =  X. 


(11) 

(12) 

(13) 


X.  =  ^ 


X 

m 


(14) 


X.  =  -i 


Y 

m 


(15) 


Xf,  =  — 
°  m 


(16) 


X,  = 


■^10 

y^sin(Xg)  + 

^  X 
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/sinlX^) 
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/2 

vz 
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The  original  version  of  the  AVCO  code  uses  a  fifth  order  Runge-Kutta  method  with  step- 
size  control  based  on  error  estimates.  However,  the  step-size  is  either  doubled  (if  the  computed 
error  is  less  than  a  specified  tolerance)  or  halved  (if  the  computed  error  is  greater  than  a  specified 
tolerance).  The  rather  arbitrary  doubling  or  halving  of  the  step-size  can  be  inefficient  because,  for 
example,  the  step-size  may  be  able  to  be  greater  than  two,  or  it  may  not  need  to  be  halved,  or  dou¬ 
bling  may  cause  the  error  to  be  too  large,  and  then  the  another  iteration  is  required,  etc.  It  is  more 
efficient  to  let  the  magnitude  of  the  computed  error  control  the  step-size.  The  code  was  modified 
to  use  a  time  stepping  method  that  computes  the  new  step-size  based  on  the  computed  truncation 
error.  The  new  method  incorporated  into  the  code  is  general  in  the  sense  that  any  order  Runge- 
Kutta  method  can  be  used.  It  is  required  only  to  input  a  table  of  coefficients  specific  to  the  order  of 
the  method.  The  Runge-Kutta  method  used  was  developed  by  Fehlberg  [8].  The  equations  for  a 
fifth  order  scheme  are  written  for  the  variable  Xj  as: 


fj.  = 

(23) 

(  "  'l 

^  x  =  o  ^ 

K  =  0,  1,2,  3,4,  5 

(24) 

K  =  0 

(25) 

7 

(26) 

K  =  0 


The  local  truncation  error  is  computed  as, 

(LTE)j  =  Xj-Xj  =  ^  =  O(dt^)  (27) 

K  =  0 

The  new  step-size  is  computed  based  on  the  local  truncation  and  is  written  as, 

1 

(  ERMAX  V 

^  ^  new  (  ^  old[^  (LTE)  )  (28) 

where  ERMAX  is  the  maximum  allowable  error  set  by  the  user. 
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Although  the  above  expressions  are  presented  for  a  fifth-order  scheme,  the  procedure  is  valid  for 
any  order  scheme.  The  user  needs  only  provide  a  set  of,  ,  and  .  Sets  of  these  coef¬ 

ficients  for  various  order  methods  are  presented  by  Fehlberg  [8]. 

3.2  Modifications  to  Newtonian  Regime  Model 

The  original  AV CO  code  uses  a  straight  Newtonian  method  to  compute  the  pressure  di  i- 
bution  on  each  of  the  faces  of  the  sabot  in  the  Newtonian  regime.  However,  straight  Newfi  .an 
does  not  account  for  the  effect  of  Mach  number,  and  tends  to  over  predict  pressure  near  stagnation 
regions.  Additionally,  there  is  an  inconsistency  in  the  code  because  the  actual  side-face  areas  of  a 
particular  sabot  are  used  to  compute  forces  and  moments  in  the  interaction  regime,  whereas  the 
side-face  areas  of  the  simple  cylindrical  shape  shown  in  Figure  2  are  used  to  compute  the  forces 
and  moments  in  the  Newtonian  regime.  This  tends  to  cause  an  overprediction  in  the  radial  dis¬ 
placement  after  the  sabot  transitions  to  the  Newtonian  regime  because  the  actual  side-face  area  for 
a  given  sabot  is  generally  less  than  that  of  the  simple  cylindrical  shape  shown  in  the  figure. 

The  code  was  changed  to  use  the  modified  Newtonian  method  for  all  Newtonian  computa¬ 
tions.  Additionally,  the  same  side-face  areas  used  to  compute  the  forces  and  moments  in  the  inter¬ 
action  region  are  now  used  to  compute  the  forces  and  moments  in  the  Newtonian  regime.  This 
results  in  a  better  prediction  of  the  radial  displacement,  especially  at  later  times,  as  will  be  shown 
in  Section  5.0. 

3.3  New  Front  Scoop  Geometry  &  Aerodynamic  Model 

The  largest  modification  to  the  modeling  was  made  to  the  front  scoop  geometry  and  asso¬ 
ciated  aerodynamic  model.  The  original  version  of  the  AVCO  code  modeled  the  front  scoop  with 
a  simple  conical  shape.  However,  experimental  designs  under  investigation  today,  for  example  the 
SLEKE  sabot,  use  a  ramp  on  the  front  scoop  for  structural  support  of  the  tip  of  the  rod.  There  is  no 
adequate  way  to  model  this  geometry  with  the  original  AVCO  code.  Additionally,  this  geometry 
alters  the  aerodynamic  flow  field  in  the  vicinity  of  the  front  scoop,  and  the  flow  field  around  this 
geometry  is  not  modeled. 

To  model  the  front  scoop  with  more  than  one  line  segment  it  was  necessary  to  develop 
force  and  moment  expressions  for  a  line  segment  arbitrarily  oriented  with  respect  to  the  point  O  in 
Figure  5.  The  basic  expressions  for  force  and  moment  are  given  by: 

/■  =  dA  (29) 

A 

ilt  =  jPhxrdA  (30) 

A 

Where  P  is  the  pressure,  h  is  the  unit  normal  vector,  and  r  is  the  moment  arm  from  a  suitably 
chosen  reference  point  (in  this  case  point  “O”).  For  the  front  scoop  h .  and  r  ■  are  given  for  line 
segment  j ,  referring  to  Figure  5  as, 
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“O” 


(d) 

Figure  5.  New  front  scoop  geometry  model. 


hj  =  COS0,  sin(t)?+ COS0 .  cos<|)}  + sin0 .  ^ 


K  =  r^sin(|)?  +  /-.cos(t)}  +  zl 


In  the  interaction  region  there  are  two  cases  to  consider.  One  is  an  attached  shock  configu¬ 
ration,  where  the  pressure  is  assumed  to  be  constant  on  the  front  scoop  behind  the  shock.  The 
other  is  a  detached  shock  configuration,  where  the  pressure  is  assumed  to  vary  linearly  in  the 
radial  direction  over  the  front  scoop. 

For  constant  pressure,  the  forces  and  moments  on  segment  j  were  derived  to  be: 


F  =  -2Psin(t)  ^)cot0 


’j±lZD 
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Finally,  it  was  necessary  to  develop  force  and  moment  expressions  for  the  front  scoop  in 
the  Newtonian  regime.  This  was  done  for  a  single  segment  in  Crimi  &  Siegelman  [2].  The  com¬ 
ponents  of  force  on  segment  j  become: 

Pxj  =  -Pco(''5+i  -'•J)  cote^  l^y^l^^sinGcos  (0(1)^-  sin <1)^ cos (j)^)  (39) 

+  |y^y,cos\.sin\] 

Py.  =  -Pco  1  -  '■;)  cote {  y^sin^e^sintt)^  +  y^y^sinG^  cosG^^tt)^  (40) 

1  2 

+  sin (j)^ cos (j)^  +  -cosGy  [y^sintj)^  (2  +  cos(^J  +  ysin(|)^]  } 
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(41) 


=  ~P«> (0+  1  - o)  +  27^y^sine^  cose^  sin(()^ 

+  yos^Qfivl  (4)^  +  sin(|)^cos(t)J  +  vjsin\]  } 

The  expression  for  the  moments  are  modified  by  an  additional  term  due  to  a  non-zero  z 
The  expressions  for  the  moments  on  segment  j  about  point  O  become: 


A^yj  =  ^0^1  -  0)  - 1  1  -  0) 

V^cose^^|- V^cos0^^sin\  +  V^sinQ _  i sin (|)^ cos 4)^11  +  F^Zj 


The  aerodynamic  model  associated  with  the  new  front  scoop  geometry  is  a  combination  of 
the  constant  pressure  and  linearly  varying  pressure  cases  described  above.  The  specific  combina¬ 
tion  depends  on  the  number  of  line  segments  used  to  model  the  front  scoop,  the  angle  of  attack  of 

the  sabot,  and  whether  the  flow  on  the  underside  of  the  sabot  is  choked  based  on  1-D  inviscid  eas 
dynamic  relations.  ^ 

For  the  two  segment  front  scoop,  two  situations  are  considered.  These  situations  are 
sketched  m  Figures  6  (a)  and  (b).  At  small  angles  of  attack,  an  attached  shock  exists  on  the  outer 
segment  and  a  detached  shock  forms  ahead  of  the  inner  segment.  The  pressure  on  the  outer  seg¬ 
ment  IS  assumed  to  be  constant  at  the  pressure  behind  the  oblique  shock.  The  pressure  on  the  inner 
segment  is  assumed  to  vary  linearly  from  stagnation  conditions  behind  a  normal  shock  corre¬ 
sponding  to  the  conditions  behind  the  oblique  shock,  to  sonic  conditions  where  the  inner  segment 
intersects  the  sabot  underside.  When  the  angle  of  attack  is  larger  than  for  an  oblique  shock  to  exist 
on  the  outer  segment,  one  detached  shock  is  assumed  to  exist  and  the  pressure  is  assumed  to  vary 
linearly  from  stagnation  at  the  outer  radius  of  the  sabot  to  sonic  at  the  inner  radius. 

For  the  three  segment  front  scoop,  the  pressure  on  each  segment  of  the  front  scoop  is 
^sumed  to  be  uniformly  at  stagnation  conditions  behind  a  normal  shock  corresponding  to  the 
incoming  freestream  conditions. 

4.0  ANALYTICAL  MODEL 

In  order  to  help  verify  the  trends  predicted  by  the  AVCO  code,  a  simple  analytical  model 
was  developed.  These  results  are  compared  to  results  predicted  by  the  AVCO  code.  A  schematic  of 
the  analytical  model  is  shown  in  Figure  7. 
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Figure  6.  Two  segment  aerodynamic  model. 


Figure  7.  Analytical  model. 


The  model  assumes  that  the  pressure  on  the  front  scoop  is  the  stagnation  pressure  behind  a 
normal  shock  corresponding  to  the  freestream  Mach  number.  This  relationship  is  well  known  and 
is  given  by, 


It  is  assumed  that  the  pressure  on  the  front  scoop  is  uniformly  at  stagnation  pressure  and  that  the 
flow  at  the  comer  point  where  the  front  scoop  intersects  the  underside  of  the  sabot  is  at  sonic  con¬ 
ditions  corresponding  to  an  isentropic  expansion  from  the  conditions  behind  the  normal  shock. 
The  flow  then  turns  through  a  Prandtl-Meyer  expansion  around  the  comer.  Since  it  is  assumed  that 
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the  conditions  at  the  comer  are  sonic,  the  Prandtl-Meyer  function  on  the  underside  is  simply 
V  (M2)  =  n-Qf.  Knowing  the  value  of  v ,  the  Mach  number  after  the  Prandtl-Meyer  expansion 
is  determined.  The  static  pressure,  ,  is  then  determined  from. 


* 

P. 


Y 


1  + 


Y-1 


Y-1 


(45) 


where  P*  is  the  pressure  corresponding  to  sonic  conditions. 

With  the  above  pressures  acting  on  the  front  scoop  and  underside  of  the  sabot,  the  forces 
and  moments  acting  on  the  sabot  petal  are  summed  assuming  that  the  sabot  pivots  about  the  aft 

end.  Summing  moments  and  applying  =  I-d^a/d  the  equation  for  rotation  is  obtained  as. 


d\t-  1 

2 

f  3  3\ 

T. 

sin  9. 

1  3  J 

(46) 


The  expression  in  the  brackets  is  the  total  moment  M^  about  the  sabot  petal  tail,  and  the  entire 
expression  on  the  right  hand  side  is  a  function  of  the  sabot  geometry  and  the  freestream  condi¬ 
tions.  The  term  F.^  is  the  same  as  given  by  equation  (33).  The  entire  right  hand  side  expression  is 

constant  for  a  given  geometry  and  freestream  conditions  so  the  equation  can  be  easily  integrated 
twice  to  obtain  an  expression  for  the  angular  velocity  and  angular  displacement  of  the  sabot  petal 
as  a  function  of  time.  The  resulting  equation  is  converted  to  a  function  of  down  range  position  by 
introducing  the  projectile  velocity  into  the  equation,  i.e.,  t  =  x/V^.  The  resulting  equation  is 
simply. 


a{x) 


Mj  x^ 


(47) 


Next,  the  forces  in  the  y  direction  are  summed  to  obtain  an  expression  for  the  radial  displacement 
of  the  sabot  as  a  function  of  down  range  position  x .  Sununing  forces  yields: 
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(P-r 

dP 


2sin(t)^ 


m 


f  2  2\ 

~  '■/ 

■ 

{r.l.  +  A^)P,-Po 

[  2  J 

cot  Gy 

COS 

2/ 

- 

O'O 


(  1  1\ 
^0-^i 


sm 


^  MpP  \ 

v^y 


{48) 


The  expression  for  angle  of  attack,  a ,  from  equation  (47)  has  been  inserted  into  (48).  This  equa¬ 
tion  can  be  integrated  twice  to  obtain  an  expression  for  the  radial  displacement  as  a  function  of 
position.  However,  to  integrate  the  equation,  the  sine  and  cosine  terms  must  be  written  in  terms  of 
their  series  expansions.  The  resulting  expression  for  radial  displacement  is  written  as: 


r{x)  = 


_  2sin(i) 


m 


(r,/.  +  A)P2-/^, 


(2  2^ 


V  2  , 


COt0. 


Ln  =  0 


(-1) 

(2n)! 


(4n+  1)  (4n  +  2) 


(49) 


-2P  (b 

O^C 


('2  2\ 

oo  „ 

Y  (-1) 

^  Mt 

2n  +  1 

V  2  j 

^  (2n+l)! 

/2  =  0 

(4n  +  3)  (4n  4-4) 

5.0  NUMERICAL  RESULTS 

The  modified  version  of  the  AVCO  code  was  used  to  compute  the  sabot  separation  trajec¬ 
tories  for  three  lAT  sabots.  The  results  are  compared  to  experimental  data  obtained  in  the  lAT’s 
Hypervelocity  Launch  Facility.  The  results  of  the  analytical  model  are  also  presented  and  com¬ 
pared  to  the  experimental  data  and  to  the  results  computed  by  the  AVCO  code.  Overall,  there  was 
good  agreement  between  the  experimental  data  and  the  results  computed  by  the  AVCO  code.  Fur¬ 
thermore,  there  was  surprisingly  good  agreement  between  the  results  computed  by  the  analytical 
model  and  the  experimental  data,  considering  the  simplicity  of  the  analytical  model. 

The  lAT  sabots  are  shown  in  Figures  8  (a)  through  (c)  and  are  referred  to  as  the  HVP_009, 
HVP_015,  and  HVP_016  sabots,  respectively.  The  HVP_009  and  HVP_015  sabots  are  similar, 
the  primary  difference  being  the  axial  location  of  the  center  of  mass.  The  front  scoop  on  these  sab¬ 
ots  is  adequately  modeled  with  two  line  segments.  The  primary  difference  between  the  HVP_016 
sabot  and  the  other  two  is  the  addition  of  a  forward  ramp  on  the  front  scoop.  This  geometry  is 
modeled  with  three  line  segments. 
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(b) 

Figure  9.  Radial  (a)  and  angular  (b)  displacement  computed  by  new  and  old  Runge-Kutta  time 
stepping  schemes. 

Before  the  modified  code  was  used  to  predict  sabot  trajectories,  the  new  Runge-Kutta  and 
time  stepping  algorithm  were  compared  to  the  old  algorithm  to  insure  that  the  new  scheme  was 
working  properly  and  producing  accurate  results. 

Figures  9  (a)  and  (b)  show  the  radial  and  angular  displacements,  respectively,  for  the 
HVP_009  sabot  computed  using  both  the  original  and  new  Runge-Kutta  and  time  stepping 
schemes.  As  previously  mentioned,  the  original  code  uses  a  5th  order  scheme.  The  results  using 
the  new  scheme  are  presented  using  2nd,  5th,  6th,  7th,  and  8th  order  schemes.  As  seen  in  the  fig¬ 
ure,  the  results  computed  by  the  new  scheme  agree  with  the  results  computed  using  the  old 
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scheme,  indicating  that  the  new  scheme  is  producing  accurate  results.  It  should  be  noted  that  the 
old  5th  order  scheme  took  189  time  steps  to  compute  the  solution  compared  to  125, 141, 152  188 
and  139  time  steps  for  the  new  2nd,  5th,  6th,  7th,  and  8th  order  schemes,  respectively,  so  the  new 
scheme  computes  the  solution  in  less  time  steps  than  the  old  one.  The  results  indicate  that  the  2nd 
order  scheme  may  be  the  most  preferable  since  it  not  only  uses  the  fewest  number  of  time  steps 
but  also  performs  the  fewest  computations  per  time  step.  ^  ’ 

Next  the  code  was  run  for  the  three  lAT  sabots  to  determine  the  overall  accuracy  of  the 
code  on  a  variety  of  sabot  geometries,  and  to  assess  the  improvements  made  to  the  original  ver¬ 
sion  of  the  code.  The  results  are  presented  and  discussed  in  the  following  paragraphs  along  with 
the  results  of  the  analytical  model. 

The  HVP_009  sabot  and  HVP_015  sabot  are  geometrically  very  similar.  The  front  scoops 
on  these  two  sabots  are  adequately  modeled  by  two  line  segments.  With  the  old  version  of  the 

code,  the  front  scoops  would  be  modeled  as  one  straight  line  segment  from  the  outer  to  the  inner 
radius. 

Both  the  modified  and  original  versions  of  the  code  were  used  to  compute  the  radial  and 
angular  separation,  and  these  results  were  compared  to  the  experimental  data  to  determine  the  rel¬ 
ative  performance  of  each  version  of  the  code.  The  experimental  data  were  taken  at  range  pres¬ 
sures  of  380  torr  for  the  HVP_009  sabot  and  at  380  torr  and  430  torr  for  the  HVP_015  sabot.  The 
computed  and  experimental  radial  and  angular  displacements  are  shown  in  Figures  10  (a)  and  (b) 
for  the  HVP_009  sabot,  and  in  Figures  1 1  (a)  and  (b)  (370  torr)  and  in  Figures  12  (a)  and  (b)  (430 
torr)  for  the  HVP_0 15  sabot.  ^ 

In  all  cases,  the  modified  version  of  the  code  computes  the  radial  displacements  more 
accurately  than  the  original  version.  Remarkably,  the  simple  analytical  model  predicts  the  radial 
displacements  more  accurately  than  either  version  of  the  code.  The  angular  displacements  com¬ 
puted  by  the  modified  and  original  version  of  the  code  are  nearly  the  same,  and  both  predict  the 
angular  displacements  accurately.  However,  the  analytical  model  does  not  predict  the  angular  dis¬ 
placements  as  well  as  the  AVCO  code. 

The  computed  and  experimental  results  for  the  HVP_016  sabot  are  shown  in  Figures  13 
(a)  and  (b).  This  sabot  was  not  modeled  using  the  original  version  of  the  code  because  there  is  no 
way  to  adequately  model  the  front  scoop  geometry  of  this  sabot  with  the  original  version  of  the 
code.  The  results  show  the  same  general  trends  as  the  HVP_009  and  HVP_015  sabots,  with  the 
computed  results  slightly  overpredicting  the  radial  displacement  while  correctly  predicting  the 
angular  displacement  and  the  analytical  model  more  accurately  predicting  the  radial  displacement 
while  somewhat  overpredicting  the  angular  displacement. 
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Figure  10.  Computed  and  experimental  radial  (a)  and  angular  (b)  displacements, 
HVP_009  sabot,  380  torn 
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(a) 


(b) 


Figure  11.  Computed  and  experimental  radial  (a)  and  angular  (b)  displacements,  sabot 
HVP_015, 370  torn 
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A 


Figure  12. 


(b) 


Computed  and  experimental  radial  (a)  and  angular  (b)  displacements,  sabot 
HVP_015, 430  torn 


Range  (in) 


(b) 

Figure  13.  Computed  and  experimental  radial  (a)  and  angular  (b)  displacements,  sabot 
HVP_016,  500  torn 


6.0  USING  THE  MODIFIED  CODE 

The  changes  in  the  input  file  for  the  modified  code  are  relatively  minor  and  relate  only  to 
the  specification  of  the  front  scoop  geometry.  Table  1  shows  the  portion  of  the  original  and  modi¬ 
fied  code  input  files  up  to  the  specification  of  the  side-face  areas.  On  line  four,  the  front  scoop 
angle  (150.0)  has  been  removed  and  replaced  with  three  additional  lines  just  before  the  specifica- 
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tion  of  the  side-face  areas.  The  first  line  contains  the  number  of  segments  used  to  d;  ibe  " 
front  scoop  geometry,  in  this  case  two.  The  second  line  gives  the  angle  of  each  segment  ing 
original  definition  of  the  front  scoop  angle.  The  third  line  gives  the  radial  location  of  t  in 
most  end  of  each  segment.  These  are  the  only  modifications  to  the  input  files  needed  t  un 
modified  version  of  the  code,  all  other  changes  are  internal  to  the  code.  The  code  (ir  .f) 
creates  an  input  file  interactively  from  user  keyboard  input  has  also  been  modified  to  ,ept 
new  front  scoop  geometry  definitions. 


Table  1.  Sample  Input,  Original  and  Modified  Code 


Original  Code  Input 


3.1251730000000002E-03  3.8773870000000003E-05  3.9285970000000001E-05 
1.4176770000000000E-06  8.8792139999999992E-07  l.OOOOOOOOOOOOOOOlE-05 
5.000000000(X)00002E-05  3.5299999999999998E-02  O.OOOOOOOOOOOOOOOOE+00 
O.OOOOOOOOOOOOOOOOE+00  2.2500000000000000E-03  1  .OOOOOOOOOOOOOOOOE-05 
45.00000000000000  137.0000000000000  150.0  0.312500000000000 

0.341670000000000 

9.8750000000000001E-03  5.7083000000000002E-02  0.15 
8530.100000000000  0.1909916700000000  -2000.000000000000 

l.OOOOOOOOOOOOOOOlE-04  1  .OOOOOOOOOOOOOOOOE-04  0.083875 
O.OOOOOOOOOOOOOOOOE+00  100000000.0000000  100000000.0000000 

100000000.0000000  100000000.0000000  l.OOOOOOOOOOOOOOOOE-04 

O.OOOOOOOOOOOOOOOOE+00  7.9000000000000006E-06  0.1083330000000000 
-0.1909916700000000  1.400000000000000  1030.150000000000 

1.1519171000000000E-03  1120.000000000000  1.000000000000000 

0  110  0  5 

0  0  1  0  0  0 
0 


Modified  Code  Input 


3.1251730000000002E-03  3.8773870000000003E-05  3.9285970000000001E-05 
1.4176770000000000E-06  8.8792 139999999992E-07  1. 0000000000000001 E-05 
5.0000000000000002E-05  3.5299999999999998E-02  O.OOOOOOOOOOOOOOOOE+00 
O.OOOOOOOOOOOOOOOOE+00  2.2500000000000000E-03  l.OOOOOOOOOOOOOOOOE-05 
45.00000000000000  137.0000000000000  0.312500000000000 

0.341670000000000 

9.875000000000000 lE-03  5.7083000000000002E-02  0.15 
8530. 1 00000000000  0. 1 9099 1 6700000000  -2000.000000000000 

1. 000000000000000 lE-04  l.OOOOOOOOOOOOOOOOE-04  0.083875 
O.OOOOOOOOOOOOOOOOE+00  100000000.0000000  100000000.0000000 

100000000.0000000  100000000.0000000  l.OOOOOOOOOOOOOOOOE-04 

O.OOOOOOOOOOOOOOOOE+00  7.9000000000000006E-06  0.1083330000000000 
-0.1909916700000000  1.400000000000000  1058.000000000000 

1.1830500000000000E-03  1120.000000000000  1.000000000000000 

0  110  0  5 

0  0  1  0  0  0 
0 
2 

90.00000000000000  160.0000000000000 

9.875000000000000  lE-03  2.8267000000000000E-02  5.7083000000000002E-02 
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7.0  CONCLUSIONS  &  RECOMMENDATIONS 


Several  modifications  were  made  to  the  AVCO  code  in  an  attempt  to  improve  its  accuracy 
mid  to  allow  the  code  to  compute  the  separation  of  more  complex  sabot  geometries.  The  two  most 
significant  modifications  were  changing  the  code  to  use  the  same  side-face  areas  to  compute  forces 
and  moments  in  the  Newtonian  regime  as  are  used  in  the  interaction  regime,  and  incorporating  a 
new  front  scoop  geometry  and  associated  aerodynamic  model  into  the  code.  The  first  modification 
results  in  an  improved  prediction  of  the  radial  displacement  over  the  original  version  of  the  code, 
whereas  the  second  allows  the  code  to  be  used  to  compute  trajectories  of  more  complex  sabot  de¬ 
signs. 

The  results  indicate  improved  predictions  over  the  original  version  of  the  code.  The  angular 
displacements  were  predicted  very  well.  However,  the  code  still  tends  to  overpredict  the  radial  dis¬ 
placements  in  all  cases  computed,  although  not  as  greatly  as  the  original  version  of  the  code.  This 
indicates  that  further  refinements  to  the  aerodynamic  model,  specifically  the  underside  pressure 
model,  may  be  needed. 

It  is  interesting  to  note  that  the  analytical  model  predicted  the  radial  displacements  very 
well  in  all  cases.  This  may  give  a  clue  as  how  to  refine  the  aerodynamic  model  in  the  code.  The 
forces  computed  with  the  analytical  model  and  the  AVCO  code  are  very  similar  except  in  one  re¬ 
spect  It  was  previously  discussed  that  the  AVCO  code  incorporates  a  “pressure  pulse”  model  into 
the  side-face  area  and  underside  force  computations  based  on  the  fact  that  previous  experiments  by 
Schmidt  [7],  on  sabots  of  the  type  shown  in  Figure  2,  indicate  that  a  pressure  pulse  occurs  on  the 
underside  and  side-face  areas  of  the  sabot  due  to  shock  reflections  and  interactions.  The  procedure 
developed  by  Crimi  &  Siegelman  [2],  approximates  the  pressure  distribution  using  an  empirical 
procedure.  The  pressure  pulse  predicted  by  this  procedure  was  compared  to  the  experimental  data 
by  Schmidt  [7],  and  shown  to  provide  acceptably  accurate  results.  However,  it  must  be  noted  that 
this  procedure  was  applied  to  sabots  with  simple  conical  front  scoops  as  shown  in  Figure  5  (a)  and 
compared  to  data  taken  using  like  sabots.  However,  it  is  not  known  what  effect  front  scoops  of  the 
type  shown  in  Figures  5  (b)  and  (c)  have  on  the  magnitude  and  distribution  of  the  pressure  pulse. 
If  the  effect  of  the  different  front  scoop  geometries  is  to  diminish  the  pressure  pulse,  then  it  is  pos¬ 
sible  that  the  code  predicts  the  angular  displacement  well  and  not  the  radial  displacement  because 
the  pressure  pulse  incorporated  into  the  code  causes  a  net  increased  lifting  force  in  the  radial  direc¬ 
tion.  But  if  the  pulse  is  centered,  or  nearly  centered,  about  the  sabot  center  of  gravity,  then  the  mo¬ 
ments  produced  by  the  forces  on  either  side  of  the  center  of  gravity  tend  to  cancel  each  other, 
thereby  having  little  or  no  effect  on  the  angular  displacement. 

In  order  to  investigate  the  effect  of  the  different  front  scoop  geometry  and  to  test  the 
hypothesis  stated  above,  an  experimental  program  on  the  sabot  shown  in  Figure  8  (c)  has  been  ini¬ 
tiated.  The  experimental  program  will  be  carried  out  in  the  Mach  5  blowdown  wind  tunnel  at  The 
University  of  Texas  at  Austin.  Detailed  pressure  measurements  on  the  front  scoop,  underside,  and 
side-face  areas  of  the  sabot  will  be  made  to  determine  whether  a  pressure  pulse  exists  for  this 
geometry  and,  if  so,  how  well  the  current  procedure  in  the  AVCO  code  predicts  its  shape  and 
magnitude.  The  data  will  be  used  to  modify  the  pressure  pulse  model  in  the  AVCO  code,  if  neces¬ 
sary.  This  will  be  documented  in  a  subsequent  report. 
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Table  2.  HVP_009  Geometric  Input  Data 
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Table  3.  HVP_015  Geometric  Input  Data 


Table  4.  HVP_016  Geometric  Input  Data 
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